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^ Abstract 
c- 1 , 

We have computed semi-inclusive spin asymmetries for proton and deuteron 
targets including next to leading order (NLO) QCD corrections and contri- 
butions coming from the target fragmentation region. These corrections have 
been estimated using NLO fragmentation functions, parton distributions and 
also a model for spin dependent fracture functions which is proposed here. 
We have found that NLO corrections are small but non-negligible in a scheme 
where gluons are polarised and that our estimate for target fragmentation 
effects does not modify significantly charged asymmetries but affects the so 
called difference asymmetries. 
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Introduction: 



Recently, the Spin Muon Collaboration (SMC) jlj] have presented a measurement of semi- 
inclusive spin asymmetries for positively and negatively charged hadrons from deep inelastic 
scattering of polarised muons on polarised protons and deuterons. This data, combined with 
previous measurements PHI] of this kind can be used to determine polarised valence and 
non-strange sea quark distributions, independently from totally inclusive data. 

Up to now, the analyses of semi-inclusive spin asymmetries have been performed 

in the naive quark-parton model, neglecting both higher order corrections and contributions 
coming from the target fragmentation region. This procedure simplifies greatly the extrac- 
tion of parton distributions and seems to be adequate given the present accuracy of the data 
and the restriction to high hadron energy fractions. 

However, taking into account that the most recent analyses of parton distributions [|]||[7|], 
which are performed in the NLO approximation from totally inclusive data, have shown the 
importance of including these effects, and that the forthcoming semi-inclusive experiments 
H promise better accuracy than the obtained so far, it is worthwhile analysing the size of 
these hitherto neglected contributions. 

Higher order corrections can be non-negligible if gluons are polarised in the proton and, 
in the case of semi-inclusive processes, require a non trivial treatment of collinear divergences 
related to the target fragmentation region. This has been addressed in references || and 
T0|| • m this last reference, the concept of fracture functions has been introduced as a mean 



to describe target fragmentation phenomena. The full NLO contributions to semi-inclusive 
cross-sections, including those related to fracture functions, have been calculated recently in 
references |]TTJ and |L_2] for unpolarised and polarised deep inelastic scattering, respectively. 

As the parton distributions and the fragmentation functions, fracture functions are non- 
pertubative objects that have to be extracted from semi-inclusive high precision experiments. 
This task is not possible yet, however one can estimate the size of the target fragmentation 
corrections effects using a sensible model for fracture functions based on parton model ideas. 
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In order to establish our notation, in the next section we show the naive quark parton 
model expressions for the semi-inclusive cross sections and the full NLO ones. We also re- 
mind the definition of the spin asymmetries in terms of the former cross sections. In the 
following section we present our choice for parton distributions, fragmentation functions and 
the main features of the model proposed for fracture functions. In the last section we show 
results and present our conclusions. 



NLO Cross Sections 



In the naive quark-parton model, the semi-inclusive cross section for the production of 
a hadron h from polarised deep inelastic scattering of charged leptons earring momentum / 



on nucleons with momentum P, is usually written as [12 



^ = AY P X>A 9jWflviW (!) 

where A is the helicity of the lepton, Cj = 4ne 2 .a 2 / x(P + I) 2 , and Aa^ denotes the difference 
between cross sections of targets with opposite helicities. This cross section is differential in 
the variables x, y and z defined by 

_ _Q^_ _P 1 q _ E h 
X ~2P-q' V ~ P-V Z ~ E N (l-x) [ ) 

where q is the transfered momentum (— q 2 = Q 2 ) and Eh and E^ are the produced hadron 
and target nucleon energies, respectively. The unpolarised cross section can easily be ob- 
tained changing the kinematical factor XYp = A(2 — y)/y for Ym = (1 + (1 — y) 2 )/2y 2 and 
removing the A's, which denote differences in polarization. Ag« is the spin- dependent parton 
distribution of flavour i and D^u is the fragmentation function of a hadron h from a parton 
i. 

It is customary to define spin asymmetries A± N , proportional to the difference between 
the number of events for antiparallel and parallel orientation of the lepton and the nucleon 
spins, which in our notation are given by 



and in the naive parton model reduce to 

A h £« e l A l( x ) D h/i(z) 
1N Ei^q(x)D h/i (z) 1 ' 

Actually, the data on these asymmetries is restricted to positively and negatively charged 
hadrons with the cross section integrated over some range of the variable z. 



The difference asymetries ||13|| , A*£ h are given by 

h+ _ h - _ Y M Aa h N + - Aa h N ~ 

A N - 777 ~h+ ~h~ K b ) 



and in this approximation have no dependence on the fragmentation functions, leading to 
expressions like 

A «+-*- = Au v + Ad v A k+-k- = AAu v ~ Ad v ^ 

for pion production on deuterium and proton targets respectively. In the next to leading 
order approximation, the polarised cross sections have the following expression 



dAa^ t—^ j f f dudp j x 



i=q,q ( 



,Z 
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A ft (-, Q 2 ) D h/l (- Q 2 ) 5(1 - u)5(l - p) 



+ A qi (-,Q 2 )D h/i (-,Q 2 )AC qq {u,p) 
u p 

+ A qi (^Q 2 )D h/g (- p ,Q 2 )AC qg (u,p) 



+ Ag(-,Q 2 )D h/i (-,Q 2 )AC gq (u,p) 
u p 



>q\U, 



+ / -{l-x) {am£(-,(1-*)z,Q 2 ) Ul-u) +AC q 
JB u L H u \ 

+ AM^,(l-x) 2 ,Q 2 )AC s (u)}} (7) 

where the AC's are the NLO coefficient functions |14[], which are proportional to a s , and 
AM/ 1 are the spin dependent fracture functions. Details about the convolution variables 
and integration limits can be found in references [11,15]. Notice that the difference between 
equations (1) and (7) is not only proportional to a s , but there is a leading order fracture 



contribution which is neglected in the most naive approximation. Obviously, the spin asym- 
metries develop much more complicated expressions, particularly, the difference asymmetries 
do not reduce just to combinations of partons distributions as in equation (6), and depend 
on the variable z. Notice that the variable z defined in equation (2) and used in equation 
(7) coincides with the one used in the analyses performed up to now, z h = P.P h /P.q || in 
the naive approximation but they differ for higher order processes, in which the hadron may 
be produced at an arbitrary angle 9 with respect to the beam direction 

1 + cosfl 

Zh = z (8) 

The z variable so defined is much more convenient for factorization purposes ||11|| . 

Inputs 

In order to feed equation (7) with parton distributions and fragmentation functions, 



we have chosen two sets of NLO parametrizations for polarised parton distributions [15 



one for unpolarised distributions [[TIJ and one for NLO fragmentation functions jH]]. The 
polarised sets reproduce the main features of the available inclusive data and are defined 



within a physical factorization presciption (MS P ), the same chosen for the coefficients in 
equation (7). In one of these sets (set 1) the gluons are polarised whereas in the other 



(set 2), the strange sea quarks are responsible for the low value of Ellis- Jaffe integral [18 
Both sets satisfy positivity constraints with respect to the unpolarised sets, something that 
is crucial for computing asymmetries. The fragmentation functions do not imply the full 
flavour symmetry relations between hadrons that were assumed in reference These 
functions were obtained as independent NLO fits to charged pion and kaon production in 
e + e~ annihilation. 

Fracture functions are a relatively new concept and have not been measured yet, so there 
are not parametrisations available for them. However, taking into account that these func- 
tions measure the probability for finding a hadron and a struck parton in a target nucleon, 



one can approximate them as a simple convolution products between known distributions. 
These are the probabilities for finding the struck parton in a nucleon carrying a fraction x of 
its momentum, the one for finding another parton in the target remnant (with momentum 
fraction constrained to the interval [0, 1 — x\) and that for its fragmentation in the observed 
hadron with momentum fraction z(l — x). Assuming that the correlation between both 
subprocesses is dominated by the momentum balance, a tipically partonic assumption, our 
proposal for the fracture functions reads as 

1 f 1 dt 

M^x,z(l-x)) = q J (x)j—- j z jJ2lk(t(l-x))D h/k (z/t) (9) 



The index j reffers to the struck parton (quark or gluon), and k denotes an intermediate 
parton which undergoes hadronization into a particle h. A sum over all possible intermediate 
flavors and momentum fractions is implied. The function N(x), given by 

N(x) = (10) 
[1 x) 

normalizes the full remnant momentum to (1— x), as required for consistency, and guarantees 
the momentum sum rule fulfilment [|l0j 

Y,JdzzM?(x,z(l-x)) = -0L (ii) 

provided 

£ [dzzD h/j (z) = l (12) 
h J 

Analogously, spin dependent fracture functions can be modelized using spin dependent 
parton distributions for the struck parton and unpolarized distributions for the remaining 
part. The normalization function is the same as in equation (10), which also guarantees the 
analogous sum rule 

Zfdzz AMj(x, z(l -x)) = (13) 
h \ L x ) 

In the next section we estimate the higher order corrections to the semi-inclusive charged 
and difference spin asymetries using the distributions presented here and our model for 
fracture functions. 



Results: 



In order to compare with the available data on semi-inclusive spin asymmetries, we com- 
pute them taking into account the production of charged pions and kaons and we integrate 
the cross sections in the variable z over the measured range. Charged kaon production adds 
negligible contributions to charged asymmetries, which are dominated by pion production, 
however we take them into account because of its role in difference asymmetries as it will 
be discussed later. 

In figure (1) we show positively (la, lb) and negatively (lc,ld) charged hadron asym- 
metries on protons using both sets for polarised parton distributions. The solid lines cor- 
respond to the most naive contribution -0(a°) and without target fragmentation effects-, 
long-dashed lines include NLO corrections to fragmentation processes, short-dashed lines 
(almost overlapping with the solid ones) takes into account fragmentation and fracture but 
at LO, finally the dotted lines (overlapping with the long dashes) show the result of the full 
NLO computation (equation 7). 

These figures show clearly that target fragmentation effects are negligible in the charged 
asymmetries for z > 0.2. This is due, at small x where target fragmentation effects are 
large, to the suppression of the full asymmetries caused by the increase of the unpolarised 
cross section. At intermediate x, the dominance of current fragmentation over target frag- 
mentation is the main reason for the smallness of the correction. At larger values of x, 
target fragmentation becomes again of the same order of current fragmentation, however 
both hadronization contributions tend to be cancelled in the asymmetry due to the fact 
that those considered here -producing spinless final states- are essentially independent of 
the initial state polarisation (that of the struck parton). The model accounts for this fact 
because it defines fracture functions in which the hadronization part is the same for the 
polarised and the unpolarised case, being the spin dependence restricted to the probability 
of finding the struck parton. In other words, 
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= (14) 

Next to leading order corrections are small but non negligible for sets with gluon polar- 
ization, as can be seen in figures (la) and (lc), if the forthcoming experiments reach the 
expected accuracy. As these corrections are dominated by those of gluon origin, they have 
no significant consequences for set 2, figures (lb) and (Id). 

The same features are observed for deuterium targets, figures (2a), (2b), (2c) and (2d). 
We also show the most recent SMC proton and deuterium data ||, and that presented by 

EMC g. 

In figures (3a) and (3b) we show the curves of figure (la) and (lc) but only for charged 
pion production against the accuracy expected from Compass || for two years running at 
100 GeV. 

It is interesting to notice that the cancellation of fracture function contributions in 
charged asymmetries allows a naive interpretation for them with less stringent cuts in z 
than those used up to now. This choice would eventually allow a substantial improvement 
of the experimental statistics. In figure (4) we show the corrections exhibited in figure (la) 
but for z > 0.1. 

A completely different situation is observed for the difference asymetries, figures (5) and 
(6), in this case calculated for values of z greater than 0.25 in order to compare with the 
experimental data presented in reference ||]. In these asymmetries the suppression due to the 
unpolarised cross sections is not present at small x so target fragmentation effects are then 
quite significant (short dashes for LO and dots for NLO) This is related to the fact that the 
asymmetries depend on the differences between the probabilities for positive and negative 
hadron production. This also affects the intermediate x region, where the differences are 
comparable in the current and target fragmentation cases. 

Corrections to equation (1) have also other serious consequences in difference asymme- 
tries. Notice that the passage from equation (5) to equation (6) implies the cancellation of 
a factor, both in the numerator and the denominator of these equations, like 



D^+ /u (z) - D w - /u (z)\ (15) 
which is found to be zero around z ~ 0.2 in different experiments [17]. The above mentioned 



corrections, however, shift the zeros of numerator and denominator in a different way causing 
large distortions (even divergencies) from the naive expectation |T3[ . These distortions make 
pointless a naive interpretation of the difference asymmetries, at least for values of z of the 
order or lower than 0.2. Notice also that in the region where 7r + production equals that of 
7T _ , kaon production acquires relevance and can not be neglected 

At variance with the charged asymmetries, difference asymmetries only allow analyses 
with a lower cut in z if large corrections are taken into account. This fact by no means 
challenges difference asymmetries. On the contrary, the comparison between the results 
comming from them and those from inclusive and the other semi-inclusive observables al- 
lows one to explore new aspects of the parton model, in particular the hypotheses of the 
model presented here. 



Conclusions: 



In this paper we have found that that NLO corrections to semi- inclusive spin asymme- 
tries, particularly those related to target fragmentation effects, are not negligible and can 
be treated quantitatively using a sensible model for fracture functions. 

Taking into account these corrections, one can safely reduce the kinematical cuts used in 
the analysis of the experimental data on charged asymmetries, correct the naive interpreta- 
tion of the difference asymmetries for z > 0.25, and give an explanation to the features of 
the data for lower z cuts. We expect that these issues will be relevant in the analyses of the 
forthcoming semi-inclusive experiments. 
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Figure Captions: 

Figure 1: Semi- inclusive asymmetries for muoproduction of charged pions and kaons on a proton 
target with z > 0.2; a) and b) for positive hadrons calculated with sets 1 and 2, respec- 
tively, c) and d) for negative hadrons. The curves correspond to the naive estimate 
(solid), adding target fragmentation effects al LO (short dashes almost superimposed 
with the previous), current fragmentation at NLO (long dashes), and the full NLO 
prediction. The data correspond to EMC and SMC experiments. 

Figure 2: The same as in figure (1) but for deuterium targets. 

Figure 3: The same as in figures (la) and (lc) but for charged pion production. The error bars 
represent the expected accuracy for two years of running of the Compass experiment. 

Figure 4: The same as in figure (la) but for z > 0.1 

Figure 5: Naive estimate of the difference asymmetry and corrections, calculated with z > 0.25 
for porton targets. 

Figure 6: The same as in figure 5 but for deuteron targets. 



11 



FIGURES 



FIGURE 1 
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FIGURE 2 




FIGURE 3 
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FIGURE 5 




